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Abstract: 

We calculate the differential decay rate for inclusive B ^ t u X transitions to order l/zn^ 
in the heavy quark expansion, for both polarized and unpolarized tau leptons. We show 
that using a systematic 1/mi, expansion significantly reduces the theoretical uncertainties 
in the calculation. We obtain for the total branching ratio BR(S -^ r u X) = 2.30 ±0.25%, 
and for the tau polarization Apoi = —0.706 ± 0.006. From the experimental measurement 
of the branching ratio at LEP, we derive the upper bound Ai < 0.8 GeV for one of the 
parameters of the heavy quark effective theory. 
December 1993 



1. Introduction 

Recently, it lias been observed that inclusive semileptonic decays of hadrons containing 
a single heavy quark allow for a systematic, QCD-based expansion in powers of l/ruQ |jl[|. 
In particular, it has been shown that the inclusive decay rates computed in the heavy quark 
limit niQ —>■ oo coincide with those obtained in the free quark decay model, while correc- 
tions of order l/mq vanish. The leading nonperturbative corrections are of order ^/toq 
and depend on only two hadronic parameters, which parameterize certain forward matrix 
elements of local dimension-five operators. These corrections have been computed for a 
number of processes 0-0- These new theoretical developments not only provide a theo- 
retical justification for the parton model, but also allow a model independent calculation 
of the nonperturbative corrections to a high level of accuracy. 

Semileptonic B meson decays into a tau lepton are potentially very interesting. First, 
they are sensitive to certain form factors which are unmeasurable with a massless lepton in 
the final state. As a consequence, the ratio of the decay rates for the tau channel and the 
light lepton channels is sensitive to the nonperturbative corrections of order l/rn^, while 
independent of the Cabbibo-Kobayashi-Maskawa matrix elements Vcb and Vub- Second, 
the possibility to study the tau polarization offers a greater variety in the experimental 
and theoretical analysis and an independent determination of parameters. Third, the 
inclusive B -^ t u X decay rate is useful for constraining or probing certain extensions of 
the standard model, such as models with many scalar fields. 

In this paper, we study the decay B ^ t u X using a combination of the operator 
product expansion and heavy quark effective theory. The main points which we address 
are the following: 

a. We calculate analytically, to order 1/m^, the total decay rate and the lepton 
spectrum for the cases of unpolarized and polarized tau leptons. Effects of the finite tau 
lepton mass are treated exactly. 

b. We study in detail the numerical predictions for the total decay rate and the tau 
polarization. In particular, we stress that by using a systematic l/mq expansion the 
theoretical predictions become more accurate than those of the free quark decay model, 
not only because l/ml corrections are included but, more importantly, because the masses 
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of the charm and bottom quarks are correlated in a specific way. 

c. We use the experimental value of the inclusive branching ratio BR(S — » r u X) to 
derive a bound on the hadronic parameter Ai, which is related to the kinetic energy of the 
6-quark inside the B meson. 

2. Analytic Expressions 

We begin by presenting the analytic results necessary for our numerical analysis. The 
techniques by which they are obtained are described in detail elsewhere 0-0, so we will 
present only a brief summary, followed by the results of the computation. 

The inclusive differential decay distribution is determined by the imaginary part of 
the time-ordered product of two flavor-changing currents, 

T^"--z [ d^xe-"'-'' {B\T {j^\x),r{Q)]\B) , (2.1) 

where J^ = C7^(l — 75)6 or J^ = 117^(1 — 75) &. Since over most of the Dalitz plot 
the energy release is large (of order mi,), the time-ordered product can be written as an 
operator product expansion, in which higher-dimension operators are suppressed by powers 
of A/?7i5, where A is a typical low energy scale of the strong interactions. To this end, 
however, it is necessary to separate the large part of the 6-quark momentum by writing 
Pb = fUbV + k, where v is the velocity of the decaying B meson. The aim is to construct 
an expansion in powers of k/mi,, where the residual momentum k is of order Aqcd- This 
separation is most conveniently performed by using the formalism of the heavy quark 
effective theory [^ . One can then evaluate the matrix elements of the resulting tower of 
nonrenormalizable operators with the help of the heavy quark symmetries. 

The leading term in the expansion reproduces the result of the free quark decay 
model , while giving a unambiguous meaning to the heavy quark mass . The leading 
nonperturbative corrections are of relative order 1/m^ and may be written in terms of two 
parameters, Ai and A2, which are related to the kinetic energy Kb of the 6-quark inside 
the B meson, and to the mass splitting between B and B* mesons [1T0(| : 



Kb^-^, m%,-ml=4X2. (2.2) 



The operator product expansion for semileptonic decays to massless leptons has been 
performed in refs. Pj-fl]. The inclusion of the tau mass in the process B ^ t u X is a 
straightforward, but cumbersome, generahzation. For the sake of brevity, we shaU present 
only the final expressions. The tau lepton can have spin up (s = +) or spin down (s = 
— ) relative to the direction of its momentum, and it is convenient to decompose the 
corresponding decay rates as 

r{B ^t{s = ±)uX) = -r±r. (2.3) 

The total rate, summed over the tau polarizations, is given by F, while the tau polarization 
is Apoi = 2F/F. The differential decay rate depends on the kinematic variables g^, Er-, and 
Ev, where q^ is the invariant mass of the lepton pair, and Ej- and E^, denote the tau and 
neutrino energies in the parent rest frame. Let us introduce a set of related dimensionless 
variables by 

q = ^, y = , X = , 2.4 

r?7,^ nib TTT'b 

and define the mass ratios 

Pj = ^ , p^ = —L^ (2.5) 

where j = c or m is a flavor label. The triple differential decay rate may be written in 
terms of five invariant form factors Wi, for which we adopt the conventions of ref. [Q. We 
obtain 



X I if - Pr)Wi + -{yx- f + pr)W2 + ^[f{y -x)-pr{y + x)]Ws (2-6) 



+ ^PTif - Pr)W^ + PrXW^ 



and 



1 -if ,. -, 2if-Pr) y. 2if-,^) _^ 



Tbdfdydx y y + ^y2_4pj \y~Vy^-'^pr J \/y^-^pr 

X I [(?' - Pr)y - 2prX]{-2Wi +W2 + X#3 + PrW^ + VW^) ^^.7) 

- (2/2 _ Apr)[xW2 + [ct - Pr)(W^ + VF5)] 



where 



W.kl^Glm^, 



1927r3 



) 



J = c or u. 



The expressions for the form factors Wi, at tree level and to order 1/m^ in the operator 
product expansion, are 



+ '^'(^){^[8'?' + 6(y + ^) - Hy + xf] + ^[8^' + 14(y + x) - 5{y + x)^ - 16]} 



"b """& 



,|,.,(,){l_5(^l_y^},,'(,,{l|(,,.),_^,,,..)-4,} 



6m? 






^^-'-f^M^^y+^^-'^^'^y^^^-''^}^'''^'^^^^^ 



w ./,^,2(Ai + 3A2) 

^4 = 5(^) —2 , 



3m^ 



(2.9) 
where 

i = l + <?^ -Pj -y-x + ie. (2.10) 

Our results for Wi, W2 and W^ coincide with those obtained in ref. Q. The two form 
factors W4 and Wr, do not contribute when the final lepton is massless, but are important 
for the tau channel. 

We now integrate over the kinematic variables q'^ and x to obtain the differential decay 
rates with respect to the rescaled lepton energy y. Notice that a subtlety specific for the 
case of non-vanishing tau lepton mass is the appearance of the second step function in 



and (|2.7| ), which gives an upper limit for the variable x. For the case of a massless 
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lepton, this limit becomes trivial (x < oo). Integrating ( |2.6| ), we find 



1 /-7T^ 

— — = 2^2/2 -4p,<j xl y" - 3y(l + p^) + 8p, 

A2 Xo 



+ Xg 



3y2 + 62/(1 + pO-12p, 



"^^(i + Pr-y) 



+ 



hxl [y^ - V(i + p^) + 22/(3 + 7pr) + 4p,(pr - 5)) 
+ 3x0 (-Sy' + 2/'(17 + Ibpr) - 2/(24 + 46p,) - p.(18p, - 70)) 

+ 3 (52/' - 102/2(1 + pr) + 4y(3 + Ap^) + 16p,(p, - 2)) 
Ai 



3m^(l + p^-2/)2 



3 (2/" - 22/3(1 + pO + 82/Pr(l + Pr) - IQpl) 
+ 2x3 (2/4 - 52/3(1 + p,) + 22/2(5 + lip, + bpl) - AQypril + p.) - 2p,(5 - 38p, + 5p2)) 
+ 3x2 ^_2^4 ^ ^y^^^ ^ p^^ _ ^2(^5 ^ 28p, + 15p2) + 522/p.(l + p^ + 18pr(l - 6p, + p^)) 

+ 6x0 (2/^ - 32/3(1 + Pr) + y2(5 + 6p, + 5p2) - 12ypr{l + p,) - 8p,(l - 4p, + p^)) 

(2.11) 



while ( |2.7|) yields 



^-y- pi 



"T OXq 



2/-2 



+ 



A2 Xo 



miil + Pr-y) 



5x2 (/ _ 2y(2 - pO + 6(1 - pO) 



+ 



+ 3x0 (-52/2 + 2/(17 - 5pO - 4(6 - 5pO) + 3 {by^ - lOy + 12(1 - p^) 
Ai 



3ml{l + pr-y)^ 



3(-2/3 + 22/2 + 42/p,-8p,) 
+ 2x3 (-2/3 + y2(5 ^ 2p,) + y(-10 - lip, + 5p2) + 6p,(5 - 3p,)) 
+ 3x^ (22/3 - 22/2(4 + 3pO + 2/(15 + 22p, - 5p2) - 24p,(2 - p,)) 



f2.12) 



+ 6x0 (-2/3 + 2/2(3 + 4p,) - 2/(5 + lip,) + 6p,(3 - p,)) 



Here 



Xo = 1 



Pj 



(2.13) 



1 + Pr - 2/ 

At this point, we like to mention that there is an elegant alternative way to obtain the 
spectra dV/dy and dV/dy. Instead of starting from a triple differential decay rate, one can 
construct the operator product expansion for the decay B -^ t + X , where the neutrino 
is now part of the final state X. The leading contribution in the expansion is given by 



a diagram containing a charm quark-neutrino loop. The discontinuity of this diagram 
gives the lepton spectrum. In this approach, the variable xq appears in a natural way as 
the upper limit in the integration over a Feynman parameter. We have checked that this 
appraoch leads to the same results as given above. 

Finally, we perform the y-integration over the kinematic region 



2y/p7 <y<l + PT-Pj 



(2.14) 



to obtain the total decay rates. For F, we find 



r 



v^Ul + 



Ai 

2ml 



1 - 7{pj + Pr) - 7(p? + pI) + p5 + p^ + p^.p,(12 - 7(p, + pj)) 



+ 



3A2 

2ml 



+ 12 1 + 



- 3 + 5(p, + pr) - l%p] + pD + 5(pf + pI) + 7p,/9,(4 - 5(p, + pj)) 
Ai + 3A2 



2ml 



^, ,^ (l+p,-l>r+^Xf + pj In (i± "^ - "' + ^>' 



4p, 



ipr 



V 2ml ' ^pTpj 



with 



(2.15) 



A = \-2{pr + pj) + {Pt-Pj] 
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(2.16) 



For r, we obtain the lengthy expression 
f 1 



6' 



[(1 - rhr) - Pj] - (1 - ^r) (3 + Ibrhr + 5rn; + m'^) 



where 



+ pj{l - mr){21 + hlrhr + 31?fi^ + llml) 

{l-rhrf 



+ p]{2l - Ibrfir - 5ml + 47m^)/(l - rhr) - 3p? 



- 2pj (3 - 3m^ - 2m;pj) In 
Ai 



Pj 



I2ml ^ [(' - ^-^' - P^ 



- (l-m^)(l + m^)^(3 + m2) 
+ pj{l + mr){2l - Grhr + 8m^ - 2ml + 11"^t)/(1 - "^r) 



+ p^(21 - 57m^ + Uml + 42ml - 99^r + 47m^)/(l - m^)^ - 3p) 



- I2p] (3 - 3m^ - 2m2pj.) In 



Pj 



M 



4m^ 



[(1 - rhrf - pj] (1 - m^)(9 + 27m^ + 70m^ + lOm^ - 15m^ - 5m^) 



— Pj(15 — 3mT- + 62m^ — 70m^ — Ahm^ — 55m^)/(l — ttIt-) 
+ pf (57 - 84m^ + d:2ml + 260?fi^ - 235m^)/(l - rhrf - ISpJ 



12pj- (8m^ - 8m^ + 3pj- + Amlpj - 15m^pj- - 10m^p|) In ^^ ^^ 

Pj 



rhr = m^r/rrib 



ipr 



(2.17) 



:2.18) 



Our results for the differential and total unpolarized decay rates confirm a very recent 



calculation of Koyrakh [|TT|. In the limit pr — > 0, these results reduce to the expressions 
given in 0-Q. In the limit Ai, A2 — > 0, corresponding to the free quark decay model, our 
results agree with ref. [l^]. In the same limit, our expressions for dV /dy and F agree with 
those of ref. |jT^. Finally, for rhr = we find —2dV/dy = dV /dy and — 2r — F, as required. 



3. Numerical Analysis 

3.1. Input parameters 

When we neglect the tiny contribution from b —>■ u transitions, the input parameters 



entering our calculations are the mass of the tau lepton [O 



m^ = 1.777 GeV, (3.1) 

the heavy quark masses m^ and rric, the hadronic parameters Ai and A2, and the quark 
mixing parameter | Vc&|. We stress, however, that not all of these parameters are indepen- 
dent. In fact, the heavy quark effective theory can be used to construct a systematic l/mg 
expansion of the masses of hadrons containing a heavy quark, in which the parameters Ai 



and A2 appear at second order. The relevant relations are []T0| 



A1+3A2 , 

rriD =mc + A h . . . , 

Zrrir 



"-c 



where we neglect higher-order power corrections. The parameter A can be associated with 
the effective mass of the light degrees of freedom in the heavy mesons 0. For each set 
of values for A, Ai, and A2, we solve ( p.2|) to find the heavy quark masses that should 
be used in the results of sect. 2. Hence, instead of the four parameters m^, tBc, Ai, A2 we 
are left with three paramters A, Ai, A2. It will turn out that the correlation between the 
heavy quark masses, which is imposed by ( |3.2| ), reduces the theoretical uncertainties in the 
results in a significant way. 

In our analysis, we use the value of A2 that is obtained from the known value of the 
B -B* mass splitting [cf. Q] 

A2 = ""'^^ ~ ""^ = 0.12 GeV^ (3.3) 

We expect this value to be accurate up to power corrections of order A/mt, ~ 10%. At this 
point, the parameters A and Ai must still be obtained from nonperturbative calculations, 
introducing a certain amount of model dependence. From a QCD sum rule analysis, one 



finds that A lies in the range |]T5| |p.6|| 

0.45< A<0.60GeV. (3.4) 



QCD sum rules have also been used to compute Ai [|T5[ [T^] [IT3] , but these calculations 



suffer from large uncertainties. There are theoretical arguments, however, that Ai should 



be negative (as one expects, since — Ai is proportional to the kinetic energy of the heavy 
quark) 0, and that its magnitude cannot be too large fl^. Here we shall use the range 



< -Ai < 0.3 GeV^ . (3.5) 

3.2. BK{B ^TuX) 

Let us then present our numerical results. It is convenient to normalize the total 
branching ratio BR(S -^ r u X) to the measured branching fraction into final states with 



an electron 20 



BR(S^ezyX) = 10.7 ±0.5 %. (3.6) 

This eliminates the otherwise significant uncertainties from the values of \Vcb\^ and m^. 
We incorporate the one-loop QCD corrections, which may be extracted from |]2I[]. Using 



as(mfe) ^ 0.22, one finds that the total rates are corrected by multiplicative factors r]j- = 
0.90 and rj^ = 0.88, respectively. What is relevant for us is the ratio, rjr/rje = 1.02, for 
which the error due to the choice of scale in the running coupling constant is very small. 
Prom ( |2.15| ) , we then obtain the main result of this section: 

BR(S^rz/X) = 2.30 ±0.25 %. (3.7) 

There are a number of points to be made regarding this result: 

a. The 1/m'^ corrections to the free quark decay model reduce the prediction for the 
rate by approximately 4%. For example, if we take the values A = 0.5 GeV and Ai = —0.25 
GeV, corresponding to the heavy quark masses rrib = 4.8 GeV and rric = 1.45 GeV, the 
spectator model result of 2.37% is lowered to 2.28%. The fractional decreases in the 
individual rates are 7.3% for the tau channel and 3.8% for the electron channel. 

b. The 11% uncertainty in our prediction ( |3.7| ) takes into account the experimental 
uncertainty in the electronic branching ratio (|3.6|) , the theoretical uncertainties in ( |3.4|) and 
( p.5| ) , as well as an estimate of higher-order power corrections which we have not included. 
The most important of these are the l/m^ corrections to the relation ( |3.2|) between mo 
and rric. We estimate the uncertainty related to this effect to be about 2%. We emphasize 
that the largest error of all of these is the experimental one, so any improvement in the 
determination of BR(S — > ei>X) will improve the accuracy of our prediction. 
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c. Numerically, the main improvement of our calculation over the spectator model 
calculation is not in the incorporation of the nonperturbative l/?n^ corrections, but in 
using rric and nib as determined by (|3.2| ) rather than treating them as uncorrelated input 
parameters. To demonstrate this, let us take for a moment the central value BR(i? -^ 
euX) = 10.7%. Then, if we allow rric and rub to vary independently within the ranges 
1.4 <m^< 1.5 GeV and 4.6 < mb < 5.0 GeV, we findB BR(S ^ r u X) = 2.27 ± 0.55%. 
However, if we use A and Ai as our input parameters and calculate rric and m^ from ( |3.2|) , 
we obtain BR(S -^ t u X) = 2.30 ± 0.09%. Accounting for the correlation between rric 
and rrib reduces the error from 24% to 4%. Even if we allow the very conservative range 
—0.5 < Ai < 0.5 GeV and 0.4 < A < 0.7 GeV (which covers a range for rric and rrib larger 
than commonly accepted), we get BR(S -^ t u X) = 2.23 ± 0.33%. This is still more 
accurate than if we allow independent variation of rric and rrib ■ 

d. We have not included the final charmless modes X^ in our numerical analysis. 
Since iKife/V'cbl ~ 0.1 [|2^, their effect on the inclusive rate for each of the electron and tau 



channel is of order 1%. The effect on the ratio of rates is even smaller and certainly well 
below our error in ( |3.7| ). 

3.3. The tau polarization 

The polarization of the tau lepton, Apoi = 2r/r, being a ratio of decay rates, is by 
itself subject to much smaller uncertainties than are the rates themselves. We find that 
the numerical value for Apoi is rather insensitive to variations of A and Ai. Allowing these 
parameters to vary within the ranges ( |3.4| ) and ( p.5| ), we find 

Apoi = -0.706 ±0.006. (3.8) 

A few points are in order regarding this result: 

a. The 1/ml corrections to the free quark decay model reduce the prediction for the 
tau polarization by approximately 4%. 



It seems to us that previous results from spectator model calculations ignore the uncertainty 
in quark masses and thus significantly underestimate the theoretical errors. 

10 



b. The 0.9% uncertainty in Apoi is due almost entirely to our estimate of higher-order 
nonperturbative corrections. The uncertainty from the theoretical errors on A and Ai is 
surprisingly small, about 0.1%. 

c. Numerically, the main improvement is again in using rric and rrih as determined by 
(P?^). As before, if we allow rric and m^ to vary independently, we find Apoi = — 0.70±0.04, 
while if we calculate them from A and Ai and allow those inputs to vary, we find Apoi = 
—0.7045 ± 0.0005. Thus, the correlation between rric and mt reduces the error on Apoi 
from 6% to 0.1%. 

d. The effect of final charmless hadronic states X^ is included in the result ( |3.8| ). The 
effect, however, is rather small: it shifts the central value of Apoi from —0.705 to —0.706, 
assuming \Vub/Vcb\ =0.1. 

e. We have also numerically calculated the tau polarization as a function of the tau 
energy, Apoi{y) = —2{dV/dy)/{dV/dy). The 1/m^ corrections to the spectator results are 
small except for a region close to the parton model endpoint ymax = 1— Pc+Pr- However, in 
this endpoint region (corresponding to |pV| ~ 1.6 GeV) our results cannot be trusted, since 
the operator product expansion becomes singular, and a resummation is necessary [^ 
12 



4. An Experimental Bound on Ai 

The only quantity concerning inclusive B ^ t u X decays which has so far been 
measured is the total branching ratio. Adding the statistical and systematic errors on the 
recent ALEPH measurement P3] in quadrature, we will use the value 



BR(S^rz/A:) = 2.76 ±0.63%. (4.1) 



Among the parameters which go into the prediction (|3.7| ), the one that is subject to the 
largest theoretical uncertainties is Ai. Using experimental data on inclusive semileptonic 
D decays (in the electron channel), one can derive the experimental lower bound [^ 



Ai > -0.5 GeV^ (4.2) 
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We find that the semileptonic B decay in the tau channel provides an upper bound on this 
quantity. Taking into account the la lower bound on the branching ratio ( f4.1| ), the la 
range for the branching ratio to electrons (|3.6| ) , our range (p.4| ) for A, and the uncertainty 
from higher-order corrections, we obtain 

Ai < 0.8 GeV^ (4.3) 



A lower bound on Ai may also be derived from the same data, but it is weaker than 
Let us clarify the role of the various input data in deriving the bound ( [4.3|) : 

a. The upper bound on Ai is sensitive to the lower bound on the total branching ratio 
in (|]1|). If, for example, BR(S ^ t u X) = 2.76% (the central value of the ALEPH mea- 
surement), then Ai < —0.7 GeV is required, which is inconsistent with (|4.2|) . Consistency 
of our theoretical analysis (which is based on the assumption of lepton universality in the 
standard model) then suggests that the branching ratio should be at the lower end of the 
range quoted by the ALEPH collaboration. This is also clear if one simply compares our 
result (|3.7| ) with the experimental result (^4.1J ). 

b. Our bound is also sensitive to the experimental value of BR(S —>■ e u X). Thus, an 
improvement in the experimental data on this mode could lead to very strong constraints 
on Ai. 

c. The bound is weakly dependent on the precise value of A. This is in contrast to 
the bound derived in [^, which is very sensitive to A through its dependence on m^. 

d. We emphasize that in deriving the bound ( [4.3| ) , we have allowed each of the input 
parameters independently to take their extreme la values, so the result we quote is rather 
conservative. 

5. Summary 

Using the operator product expansion and heavy quark effective theory, we have cal- 
culated the inclusive rate for B meson decays into a tau lepton plus anything, including 
the effects of the non-vanishing tau mass and the leading nonperturbative corrections of 
order 1/ml. For the total branching ratio, we find BR(S —>■ t i> X) = 2.30 ± 0.25%. For 
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the polarization of the tau lepton, we obtain Apoi = —0.706 ± 0.006. The effect of the 
l/m^ corrections on both observables is a smaU shift of about 4% as compared to the free 
quark decay model. However, we have emphasized that, from the numerical point of view, 
the most significant improvement of our analysis is to implement the tight correlation be- 
tween the heavy quark masses rub and rric, which is imposed by the structure of the heavy 
quark expansion. This correlation reduces the theoretical uncertainties in a very significant 
way. The experimental value for the branching ratio allows us to derive the upper bound 
Ai < 0.8 GeV on one of the fundamental parameters of the heavy quark effective theory. 
While this paper was in writing, we received a preprint by Balk et al. [^ with results 



similar to ours and to ref. |]Tl|. Our study of the tau polarization in sect. 2, and the 



numerical analysis of sects. 3 and 4, have no overlap with either of these papers. 
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